Appln. No. 10/813,563 

Preliminary Amendment dated August 12, 2004 
Reply to Notice to File Missing Parts of June 15, 2004 

This listing of claims will replace all prior versions, and listings, of claims in the application: 
Listing of Claims: 

Claim 1 (original): A method for producing three-dimensional information of an object (4) in 
medical X-ray imaging, characterized in that 

the object is X-radiated from at least two different directions and the 
said X-radiation is detected to form projection data of the object (4) 

the object is modelled mathematically utilizing the projection data to 
solve the imaging geometry and/or the motion of the object, where the solving concerns 
either some or all parts of the imaging geometry and/or the motion of the object. 

and said projection data and said mathematical modelling of the object 
are utilized in Bayesian inversion based on Bayes' formula 

p(m) 

to produce three-dimensional information of the object, the prior distribution ppr(0) 
representing the prior knowledge of the imaging geometry and/or the motion of the object 
(4), the prior distribution ppr(x) representing mathematical modelling of the object, x 
representing the object image vector, which comprises values of the X-ray attenuation 
coefficient inside the object, 9 representing the parameter vector of the imaging geometry 
and/or the motion of the object (4), m representing projection data, the likelihood distribution 
p(m|x,0 ) representing the X-radiation attenuation model between the object image vector x, 
geometry parameter vector 0, and projection data m, p(m) being a normalization constant and 
the posteriori distribution p(x,0|m) representing the three-dimensional information of the 
object (4) and the imaging geometry including the motion of the object. 
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Claim 2 (original): A method according to claim 1, characterized in that the three- 
dimensional information of the object (4) is one or more two-dimensional images 
representing X-ray attenuation coefficient along slices through the object. 

Claim 3 (original): A method according to claim 1, characterized in that the three- 
dimensional information of the object (4) is a three-dimensional voxel representation of the 
X-ray attenuation in the object. 

Claim 4 (original): A method according to claim 1, characterized in that the measurement 
model is m = A0x+e, where matrix AG contains the weights describing how much each voxel 
contributes to the X-ray attenuation along the X-ray paths and the noise e is independent of 
object image vector x and the geometry parameter vector 9 leading to the likelihood 
distribution 

p(m \x,6) = p noise (m - A e x) . 

Claim 5 (original): A method according to claim 1, characterized in that the said 
mathematical modelling comprises that X-radiation attenuates when passing the object (4), 
which means that every image voxel is nonnegative. 

Claim 6 (original): A method according to claim 1, characterized in that mathematical 
modelling is expressed by the formula: 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the value U N (x) 
depends only on the values of voxels belonging to the neighborhood N, and a is a positive 
regularization parameter used to tune the width of the prior distribution. 
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Claim 7 (original): A method according to claim 1, characterized in that the 3D tomographic 
problem is divided into a stack of 2D tomographic problems and on every such 2D problem, 
the mathematical modelling is expressed by the formula: 



where the sum is taken over a collection of 2D neighbourhoods N and the value U N (x) 
depends only on the values of pixels belonging to the neighborhood N, and a is a positive 
regularization parameter used to tune the width of the prior distribution, and the 2D 
tomographic problems are related to each other by the formula 



where the sums are taken over all pixels (k=l,...,K, q=l,...,Q ) and y > 0 is another 
regularization parameter. 

Claim 8 (original): A method according to claim 7, characterized in that the neighborhoods 
consist of two adjacent pixels and U calculates a power of the absolute value of the 
difference, leading to the formula 



where s is a positive real number. 

Claim 9 (original): A method according to claim 8, characterized in that s=l corresponding 
to total variation (TV) prior describing objects (4) consisting of different regions with well- 
defined boundaries. 

Claim 10 (original): A method according to claim 1, characterized in that mathematical 
modelling is qualitative structural information of the target where the structural information 



p pr (x) =exp^aYy N (*)) 



N 



P r3D(xG)) = exp(- Y EZ|xG)[k,q]-xG-l)[k,^ |), 



f K-\ Q K Q-\ 
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is encoded in prior distributions that are concentrated around object image vectors x that 
correspond to the physiological structures of the object (4). 

Claim 11 (original): A method according to claim 1, characterized in that mathematical 
modelling consists of a list or probability distribution of possible attenuation coefficient 
values in the object (4). 

Claim 12 (original): A method according to claim 1, characterized in that the X-ray imaging 
geometry, such as X-ray source position, has unknown error modelled in the distribution 
p(m|x,9). 

Claim 13 (original): A method according to claim 1, characterized in that the X-radiation 
measurement noise is Poisson distributed. 

Claim 14 (original): A method according to claim 1, characterized in that the medical X-ray 
imaging is dental radiography. 

Claim 15 (original): A method according to claim 1, characterized in that the medical X-ray 
imaging is surgical C-arm imaging. 

Claim 16 (original): A method according to claim 1, characterized in that the medical X-ray 
imaging is mammography. 

Claim 17 (original): A method according to claim 1, characterized in that three-dimensional 
information of the object (4) is produced on the basis of the maximum a posteriori estimator 
(MAP) which is calculated by the equation: 

P(y map I m ) = maxpO | m) 



Page 9 of 18 



Appln. No. 10/813,563 

Preliminary Amendment dated August 12, 2004 
Reply to Notice to File Missing Parts of June 15, 2004 

m representing projection data and y=(x,0) representing the vector composed of the object 
image vector and the geometry parameter vector and where the maximum on the right hand 
side of the equation is taken over all y. 

Claim 18 (original): A method according to claim 1, characterized in that three-dimensional 
information of the object (4) is produced on the basis of the conditional mean estimator 
(CM), which is calculated by the equation: 

y C M = jyp(y I m)dy 

where m represents projection data and y=(x,8) represents the vector composed of the object 
image vector and the geometry parameter vector. 

Claim 19 (currently amended): A medical X-ray device (5) arrangement for producing three- 
dimensional information of an object (4) in a medical X-ray imaging, characterized in that 
the medical X-ray device (5) arrangement comprises: 

an X-ray source (2) for X-radiating the object from at least two 

different directions 

a detector (6) for detecting the X-radiation to form projection data of 

the object (4) 

means (15) for modelling the object (4) mathematically utilizing the 
projection data to solve the imaging geometry and/or the motion of the object, where the 
solving concerns either some or all parts of the imaging geometry and/or the motion of the 
objectr 

and the medical X-ray device (5) arrangement includes compri s e s 
means (15) for utilizing said projection data and said mathematical modelling of the object in 
Bayesian inversion based on Bayes' formula 
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p{x,0 | m) = 



_P P r(0)P P r(x)P(™\x>0) 



p(m) 



to produce three-dimensional information of the object, the prior distribution p pr (0) 
representing the prior knowledge of the imaging geometry and/or the motion of the object 
(4), the prior distribution p pr (x) representing mathematical modelling of the object, x 
representing the object image vector, which comprises values of the X-ray attenuation 
coefficient inside the object, 0 representing the parameter vector of the imaging geometry 
and/or the motion of the object (4), m representing projection data, the likelihood distribution 
p(m|x,0 ) representing the X-radiation attenuation model between the object image vector x, 
geometry parameter vector 0, and projection data m, p(m) being a normalization constant and 
the posteriori distribution p(x,0|m) representing the three-dimensional information of the 
object (4) and the imaging geometry including the motion of the object. 

Claim 20 (original): A medical x-ray device (5) arrangement according to claim 19, 
characterized in that the three-dimensional information of the object (4) is one or more two- 
dimensional images representing X-ray attenuation coefficient along slices through the 
object. 

Claim 21 (original): A medical x-ray device (5) arrangement according to claim 19, 
characterized in that the three-dimensional information of the object (4) is a three- 
dimensional voxel representation of the X-ray attenuation in the object. 

Claim 22 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the measurement as 

m = A 0 x+e, 
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where matrix A e contains the weights describing how much each voxel contributes to the X- 
ray attenuation along the X-ray paths and the noise e is independent of object image vector x 
and the geometry parameter vector 0 leading to the likelihood distribution 

p( m \x,0) = p noise (m - A G x) . 

Claim 23 (original): A medical X-ray device (5) arrangement according to claim 19 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically so that X-radiation attenuates when passing the 
object (4), which means that every image voxel is nonnegative. 

Claim 24 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically by the formula: 

P pr to =exp(na2 c/ * <») 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the value U N (x) 
depends only on the values of voxels belonging to the neighborhood N, and a is a positive 
regularization parameter used to tune the width of the prior distribution. 

Claim 25 (original): A medical x-ray device (5) arrangement according to claim 19, 
characterized in that the 3D tomographic problem is divided into a stack of 2D tomographic 
problems and on every such 2D problem, and the medical X-ray device arrangement 
comprises means (15) for modelling the object (4) mathematically by the formula: 
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P pr (*) =exp£a2Y/* (*)) 



N 



where the sum is taken over a collection of 2D neighbourhoods N and the value U N (x) 
depends only on the values of pixels belonging to the neighborhood N, and a is a positive 
regularization parameter used to tune the width of the prior distribution, and the 2D 
tomographic problems are related to each other by the formula 

pr3D(x(j)) = exp(- y EE | xG)[k,q] - x(j-l)[k,q] |), 

where the sums are taken over all pixels (k=l,...,K, q=l,...,Q ) and y > 0 is another 
regularization parameter. 

Claim 26 (currently amended): A medical X-ray device (5) arrangement according to claim 
25, characterized in that the neighborhood systems comprise con s i s t of t wo neighboring 

pixels xj, xk or voxels xj, xk and ^ N ^ calculates a power of the 



( K-\ Q K Q-l 

^ k=\ q=l k=l g=l 



absolute value of the difference, leading to the formula where s is a positive real number and 
a is a regularization parameter used to tune the width of the prior distribution. 

Claim 27 (original): A medical X-ray device (5) arrangement according to claim 26, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically by setting s=l corresponding to total variation (TV) 
prior describing objects consisting of different regions with well-defined boundaries. 
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Claim 28 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically by assuming that mathematical modelling is 
qualitative structural information of the target where the structural information is encoded in 
prior distributions that are concentrated around image vectors x that correspond to the 
physiological structures of the target. 

Claim 29 (currently amended): A medical X-ray device (5) arrangement according to claim 
19, characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically by assuming that mathematical modelling comprises 
con s i s t s of a list of possible attenuation coefficient values in the object. 

Claim 30 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically by assuming that the X-ray imaging geometry, such 
as X-ray source position, has unknown error modelled in the distribution p(m|x,0). 

Claim 31 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
modelling the object (4) mathematically by assuming that X-radiation measurement noise is 
Poisson distributed. 

Claim 32 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is dental radiography. 

Claim 33 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is surgical C-arm imaging. 
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Claim 34 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray imaging is mammography. 

Claim 35 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
producing three-dimensional information of the object (4) on the basis of the maximum a 
posteriori estimator (MAP), which is calculated by the equation: 



m representing projection data and y=(x,9) representing the vector composed of the object 
image vector and the geometry parameter vector, and where the maximum on the right hand 
side of the equation is taken over all y. 

Claim 36 (original): A medical X-ray device (5) arrangement according to claim 19, 
characterized in that the medical X-ray device arrangement comprises means (15) for 
producing three-dimensional information of the object (4) on the basis of the conditional 
mean estimator (CM), which is calculated by the equation 



where m represents projection data and y=(x,9) represents the vector composed of the object 
image vector and the geometry parameter vector. 



/?Omap I m ) = max/?0 1 m) 
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